function ODcalc = ODcalculate2D(fitparams,coormatrix,lreshape)

    ODcalcvec=gauss2derror(fitparams,coormatrix);    %OD calculated from fit paramters
    ODcalc = reshape(ODcalcvec, size(lreshape));
end